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Selected recent results from experimental studies of radiative penguin decays of B mesons by the Belle and 
BABAR collaborations are discussed: preliminary findings from a first inclusive measurement of b — > sj using a 
hadronic tag by BABAR, first presented at this conference, updated preliminary results Belle on B pAy, first 
shown at the Moriond QCD workshop earlier this year, and a recently published BABAR study of B — > and 
B ljj decays. 
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1. Introduction 

Within the standard model of particle physics 
(SM), the rare, flavor-changing- neutral-current decays 
b — > c?7 and b — > sj are forbidden at tree level. The 
leading-order processes are one-loop electroweak pen- 
guin diagrams as shown in Figure 1, where the top 
quark is the dominant virtual quark contribution. In 
the context of theories beyond the SM, new virtual 
particles may appear in the loop, which could lead to 
measurable effects on experimental observables such 
as branching fractions and CP asymmetries [1]. 
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Figure 1: Example leading-order Feynman diagram for 
b — > sj and b — > dj transitions. 



The shape of the photon energy spectrum is insen- 
sitive to non-SM physics [2] but can be used to de- 
termine the Heavy Quark Expansion (HQE) param- 
eters Mssand /i^, which are related to the mass and 
momentum of the b quark within the B meson. Im- 
proved measurements of these parameters are relevant 
to, e.g., reduce the error in the CKM matrix ele- 
ments \V c b\ and \V u b\ determined from semi-leptonic 
5-meson decays [3]. 

In the following, three recent experimental studies 
of radiative penguin decays are discussed 



2. Measuring b — ► 57 with a hadronic tag 

Previous inclusive measurements of b — » 87 transi- 
tions, by the CLEO [4], Belle [5] and BaBar [6, 7] 
collaborations, used either fully inclusive event sam- 
ples, requiring only a high-energy lepton in the event 
or no tag at all, or a combination of several exclusive 
decay modes. A new (preliminary) analysis, based on 
BABAR data corresponding to an integrated luminos- 
ity of 210 fb _1 , uses for the first time a recoil method 



to select candidate B decays, where BB events are 
tagged by a fully reconstructed hadronic decay of one 
the B mesons, in the following referred to as the tag 
B, and radiatively decaying signal B mesons are re- 
constructed from the remaining particles in the event. 
While the hadron tag efficiency is low (about 0.3%), 
this method allows for the inclusive study of b — > 57 
decays in a relatively clean environment and a deter- 
mination of the momentum, charge, and flavor of the 
B mesons. Thus, it is possible to measure the photon 
spectrum in the rest frame of the signal B, to separate 
charged and neutral B mesons and to determine the 
CP-asymmetry Ac p. 

In order to reconstruct a large sample of tag B 
mesons, hadronic decays to DY ± and D*Y ± final 
states are selected, where Y ± denotes relevant combi- 
nations of mesons - 7r ± , K ± , Kg, and 7r° - with a total 
charge of ±1 [8]. Those particles in the event that are 
not reconstructed as part of the tag B are required 
to include an isolated, well-reconstructed photon can- 
didate with an energy E^ > 1.3 GeV in the signal B 
meson rest frame. Photon candidates that are found 
to be consistent with stemming from 7r° or 77 decays 
are rejected. Continuum backgrounds (e + e~ — > qq, 
with q = u,d,s,c) are suppressed using a Fisher dis- 
criminant combining twelve quantities that are sensi- 
tive to event shape differences between B decays and 
continuum processes. 

The numbers of remaining B and non-5 events 
are determined by means of fits to the beam-energy- 
substituted mass Mes — V s /^ ~ Pb- 1 The hand 
side of Figure 2 shows the preliminary results as a 
function of the photon energy. The points are from 
the data; the solid histogram was obtained from a BB 
Monte Carlo sample (excluding the signal decay B — > 
X s j) and then scaled according to the results of a fit 
to the data in the region 1.3 GeV < E 1 < 1.9 GeV. 
For E 1 > 1.9 GeV 119 ±22 B -> X s j signal events are 
observed over a BB background of 145 ± 9 events. 

The differential decay rate (1/Y b)(o]T / dE^ is mea- 



-'^Here, y/s is the total energy in the center of mass (CM) 
frame, and ps denotes the B candidate CM momentum. 
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Figure 2: Left: the measured numbers of B events as a function of photon energy. The points are from data; the 
histogram is from a BB Monte Carlo sample which excludes the signal decay B — > X s j; Right: the partial branching 
fractions (1 /T B)(dT / dE^ with statistical and total error. All results are preliminary 



Table I Preliminary results for the b — > sj branching fraction and moments of the photon energy spectrum, with 
statistical and systematic errors, for different minimum photon energies ii/cut • 



^cut(GeV) 


B(B X sl ) (10- 6 ) 




(£ 7 >(GeV) 


((#7 


-<£ 7 )) 2 )(GeV 2 ) 


1.9 


366 ± 85 


zb 59 


2.289 


zb 0.058 zb 0.026 


0.0334 


zb 0.0124 zb 0.0263 


2.0 


339 zb 64 


zb 47 


2.315 


zb 0.036 zb 0.020 


0.0265 


zb 0.0057 zb 0.0203 


2.1 


278 zb 48 


± 34 


2.371 


zb 0.025 ± 0.011 


0.0142 


=b 0.0037 =b 0.0111 


2.2 


248 zb 38 


zb 26 


2.398 


zb 0.016 =b 0.006 


0.0092 


=b 0.0015 =b 0.0061 


2.3 


207 zb 30 


zb 19 


2.427 zb 0.010 zb 0.007 


0.0059 


zb 0.0007 =b 0.0073 



sured in bins of the signal B rest frame photon energy 
in the range 1.9 < E 1 < 2.6 GeV. For a given bin i 
this is determined according to 

T B dE 1 SiN B C tag ' 1 j 

where N{ is the number of B events in the bin, b{ 
is the number of B mesons from decays other than 
b — > (s, d)7, and Nb is the total number of B mesons 
in the sample. The efficiency Si corrects for both ac- 
ceptance and bin-to-bin resolution effects, and the cor- 
rection factor Ctag accounts for any dependence of the 
hadronic B tag probability on the the presence of a 
B — > X s 7 final state. 

In the right part of Figure 2 the partial branching 
fraction is shown after all corrections. The systematic 
uncertainties - mainly arising from the modelling of 
the BB background, the Mes fit parametrization, the 
description of the detector response, and the depen- 
dence on the B — > X s j signal model - are included 
here. The results for the integrated branching fraction 
and moments of the photon energy spectrum above 
different minimum photon energies i^ cu t are summa- 



rized in Table I. Figure 3 shows the first and second 
moments of the photon energy spectrum as a function 
of £?cut; good agreement with previous measurements 
is found. 

Using [3] to extrapolate the measured branching 
fraction down to a minimum photon energy of 1.6 GeV 
yields 

B(b^ s 7 )k>i.6GeV = (3.9=b0.9=b0.6) x lO" 4 . (2) 

This is in good agreement with the current experi- 
mental world average [12] of (3.55 =b 0.26) x 10 -4 as 
well as recent NNLO QCD calculations [9-11]. 

All these preliminary results are currently limited 
by statistics; adding more data in future measure- 
ments is expected to significantly reduce also the sys- 
tematic uncertainties. It should also be noted that 
the recoil method is complementary to those of other 
measurements of b — > 57 transitions; the largely inde- 
pendent systematic uncertainties will facilitate a com- 
bination of the results. 
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Figure 3: First and second moments of the photon energy spectrum as a function of E C ut> For comparison, the results 
from [6] and [7] are show as well. 



3. Study of B pA^y decays 

The decay B + — > pAj was first observed by Belle 
in 2005 based on a data sample corresponding to an 
integrated luminosity of 140 fb _1 [13]. These results 
have now been updated with a preliminary measure- 
ment 2 using 414 fb _1 of Belle data. This is a rare 
process; the SM branching fraction is predicted to be 
« lx 10 -6 . The final state baryons restrict the phase 
space, thus giving access to low photon momenta. The 
experimental study of this decay provides information 
on the baryon production mechanism. For example, 
the spin of the s quark can be probed through a helic- 
ity analysis; large wrong- helicity contributions would 
indicate physics beyond the SM. 

Figure 4 shows the projections of a two-dimensional 




0.2 0.4 

AE (GeV) 



m ES (GeV) 



Figure 4: The AE (left) and Mes (right) distributions 
for selected B + — > pA^y in the pA invariant mass region 
< 2.8 GeV, together with the result of the fit (solid curve), 
which is also shown separated into signal (dash-dotted) 
and background (dashed) components. 



2 Since this conference, this analysis has been finalized and 
published [16]. 



fit to the distributions of the energy difference AE = 
E* B - ^beam' wnere is tne CM energy of the B me- 
son candidate and EJ eam is the CM beam energy, and 
Mes for selected events with an invariant di-baryon 
mass m p A < 2.8 GeV. About 100 signal events are 
found, corresponding to a statistical significance of 
14<j. Turning this into a branching fration and ex- 
trapolating to the full di-baryon mass range yields 



B = (2.45±g;|| ± 0.22) x 10" 



(3) 



somewhat above current SM predictions [14, 15]. 

A near-threshold enhancement is found in the di- 
baryon mass distribution, see the left part of Figure 5. 
This was previously observed for other baryonic B de- 
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Figure 5: Measured distributions of the pA invariant mass 
(left), compared to theoretical expectations, and cos(# p ), 
where the helicity angle P is the angle between proton 
and photon direction in the di-baryon rest frame. 

cays and has since then been the focus of considerable 
theoretical interest. 

The right hand side of Figure 5 shows the mea- 
sured distribution of cos(# p ), where the helicity an- 
gle P is the angle between proton and photon di- 
rection in the di-baryon rest frame. It is consistent 
with a short-distance b — > sj description of this de- 
cay. The asymmetry with respect to cos(0 p ) = is 
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Table II Measured branching fractions B, angular asym- 
metries Ao, and CP asymmetries Acp- All results are 
preliminary. 



Mode 



B + -> pA-y 
B + -h. pAn 
B° -> Mtt- 



5 (io- 6 ) 



2.45l°;3g =b 0.22 +0.29 ±0.14 +0.17 ±0.17 
3.00lo.53 ± 0.33 -0.16 ±0.18 ±0.01 ±0.17 
3-23±g;|| ± 0.29 -0.41 ±0.11 -0.02 ±0.10 



determined to be A e = 0.29 ±0.14. The CP asymme- 
try, Acp = 0.17 ±0.17, is found to be consistent with 
zero. 

Table II summarizes these results and includes, for 
comparison, corresponding measurements for the de- 
cays B + — > pAir , which is observed for the first time, 
and B° -> pAn~ . 



4. B — > p/^7 branching fractions 

With respect to 6 — > 87, the corresponding 6 — > ^7 
branching fractions are suppressed by a factor w 0.04; 
the SM branching fractions for the decays of 5 mesons 
to p7 and coj are (9(10 -6 ). While the calculations of 
the exclusive decay rates have large uncertainties due 
to non-perturbative long-distance QCD effects, some 
of this uncertainty cancels in the ratio of B — > p/ooj to 
B — » if* 7 branching fractions. Since the dominant di- 
agram involves a virtual top quark, this ratio is related 
to the ratio of Cabibbo-Kobayashi-Maskawa (CKM) 
matrix elements \V td /V ts \ [17, 19, 28] via 3 
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Figure 6: AE and Mes projections of the fits for the decay 
modes B + -> p + <y (top), B° p°j (middle), and B° 
ujj (bottom). In each plot, the signal fraction is enhanced 
by selections on the other fit variables. The points are 
data, the solid line is the total of all contributions, and 
the long- dashed (dashed- dotted) line is background-only 
(signal-only) . 



B[B -> (p/coh] 



V td 




[1+AR]. 

(4) 



Physics beyond the Standard Model could affect these 
decays, creating inconsistencies between this measure- 
ment of |Vtd/Vts| and that obtained from the ratio of 
B° and B® mixing frequencies [20]. 

After previous searches by BABAR [21] and CLEO 
[22], which found no evidence for the decays B — > pj 
and B — > CJ7, an observation of the decay B° — > p°j 
was reported by the Belle collaboration [23]. A re- 
cently published BABAR study [24] of the decays B + — > 
p + 7, B° — > p°7, and B° — > cc>7 is summarized in the 
following. 

From a data sample containing 347 million BB 
pairs, which corresponds to an integrated luminos- 
ity of 316 fb _1 , the decays B — > pj and B — > cwy 
are reconstructed by combining a high-energy photon 



3 The coefficient £ is the ratio of the form factors for the 
decays B — > and B — ► K*j and AR accounts for different 
dynamics in the decay. 



with a vector meson reconstructed in the decay modes 
p° — > 7T + 7r _ , p + — > 7T + 7r°, and a; — > 7r + 7r _ 7r°. The 
dominant source of background is continuum events 
(e + e~ — > qq, with # = d, s,c) that contain a high- 
energy photon from 7r° or 77 decays. Other back- 
grounds include photons from initial-state radiation 
processes, decays of B — > if*7 (if* — > if7r), decays 
of B — > (p/u)tt or 5 ^ (p/uo)rj and combinatorial 
background from higher- multiplicity b — > 87 decays. 

High-energy photons 7 for which the combination 
with another photon 7' is found - based on a like- 
lihood ratio constructed from E 1 > and ra 77 / - to be 
consistent with a 7r° or 77 decay are rejected. A neural 
network, exploiting differences between background 
and B decays in, e.g., lepton and kaon production 
(through flavor-tagging variables described in [25]) 
and event shape, is used to suppress continuum back- 
ground events. 

The signal content of the data is determined by a 
multidimensional unbinned maximum likelihood fit, 
which is constructed individually for each of the three 
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Table III The signal yield (n s i g ), significance (U) in standard deviations including systematic errors, efficiency (e), and 
branching fraction (£>) for each mode. The errors on n s i g are statistical only, while for the branching fraction the first 
error is statistical and the second systematic. 
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signal decay modes. All fits use AE, Mes, cos Oh, 4 
and the neural network output N. For decays B° — > 
CJ7 (a; — > 7T + 7r - 7r°), the cosine of the angle between 
the 7r + and 7r° momenta in the 7r + 7r _ rest frame 
(Dalitz angle) is added as a fifth observable. In the fit, 
signal, continuum background, B — > K*^ decays, and 
other B backgrounds are considered as hypotheses for 
the origin of the events. 

Figure 6 shows the data points and the projections 
of the fit results for AE and Mes separately for each 
decay mode. The signal yields and corresponding 
branching fractions are listed in Table III. 

The isospin symmetry is tested by measuring the 
quantity T(£+ -> p+ 7 )/[2r(£° -> p° 7 )] - 1 = 
—0.35 ± 0.27, which - within the large uncertain- 
ties - agrees with the theoretical expectation [17]. 
The measured branching fractions are consistent with 
the isospin relation Y B ^ p + 1 = 2T B ^ p o 1 = 2Tb^uj 
among the widths of the individual decays. Impos- 
ing this relation as a constraint, the isospin- averaged 
branching fraction is determined from a simultaneous 
fit to the three decay modes to be 

B[B - (p/wft] = (1.25±g:it ± 0.09) x 10" 6 ; (5) 

the significance of the signal is 6.4<j, including system- 
atic uncertainties. Excluding the B° — » ivy mode from 
the simultaneous fit yields B(B — > pj) = (1-36^27 ^ 
0.10) x 10- 6 . 

As shown in Table III and Figure 7, all these results 
agree well with corresponding measurements by the 
Belle collaboration, published in [23]. 5 
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Figure 7: Measured B — > pj and B — > uoj branching frac- 
tion, compared to theoretical predictions. 



Using the world average value of B(B — > 
Equation (4), and theory input from [27], this trans- 
lates into 

\V td /V ts \= 0.200±° f 2 l±0Mb, (6) 

where the first error is experimental and the second 
is theoretical. This result is in very good agreement 
with the measurement of this ratio from the study of 
B° and B° mixing [20]. 
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4 The helicity angle Oh is denned as the angle between the 
B momentum vector and the tt~ track calculated in the p rest 
frame in the case of a p meson, or the angle between the B 
momentum vector and the normal to the uj decay plane for an 
u) meson. 

5 It should be noted that since this conference Belle has re- 
ported updated results based on a significantly larger dataset 
(« 600 fb -1 ). Good agreement with the previous measure- 
ments is found; still no significant signal is seen for the decay 
B° — > ojj. Details can be found in [26]. 
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